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Abstract 

The statistics of random band-matrices with width and strength of the 
band slowly varying along the diagonal is considered. The Dyson equation for 
the averaged Green function close to the edge of spectrum is reduced to the 
Painleve I equation. The analytical properties of the Green function allow to 
fix the solution of this equation. The former appears to be the same as that 
arose within the random-matrix regularization of 2d-gravity. 

1. The random band-matrices are now usually considered as the theory of quasi 
ID systems. For example it is thought that band matrices may adequately depict the 
properties of electrons in thin wires (see also p], |3| ) . 

In this note we would like to consider the natural extension of the usual band- 
matrix model. We will investigate the statistical features of quasi-banded matrices 
with the parameters of the band (width, strength of interaction etc.) slowly varying 
along the diagonal. 

For application to condensed matter physics or energy level statistics in complex 
quantum systems it may be also useful to add some regular diagonal part to the 
random Hamiltonian (see [|J). Thus we deal with the Hamiltonian 

H = H° + V , H% = fc(i)<5 y , (1) 

where h(i) is some smooth monotonic function and the elements of random matrix 
Vij are Gaussian distributed. 

As usual, due to the Whick theorem, the Gaussian ensemble may be defined by 
the second moment: ^ _ 

VijVrnn J i&jm&in ~\~ ^jn^im) • (2) 

Here F(i,j) = F(j, i) and b ^> 1 is the typical width of the band. We are interested 
in the quasi band case. Therefore we suppose that F(i,j) decreases rather rapidly 



when \i — j\ > b, but is a very slow smooth function of % + j. In particular the sum 
over j 

Wf{i )] ) = W{t)~1 (3) 

j 

is a very smooth function of i, as well as the regular function h(i), 

dW dh 1 1 
6— ~ b— ~ - < 1 . (4) 

di di d 

The correlator (FJ) describes the Gaussian ensemble of real symmetric matrices. 
One can also consider the ensemble of general Hermitean quasi band-matrices. To 
this end it is enough to withdraw the second term in brackets in (0). All the results 
obtained in this paper are relevant for both ensembles. 

2. The "physical quantity" we would like to consider is the Green's function 
G = {E-H)- 1 . As follows from (g) 

G~ = G(i)5 l3 . (5) 

The formal expansion in a series in V/E allows to write down the Dyson-type equa- 
tion: 



^ij + X/ — ~~i — ~(GkkGij + GkiGkj) 



E-h(i 



Here we have performed Whick contraction of the first element in each (V) n with 
any other and then resumed the series. In order to simplify this exact equation it 
is easy to note that the second term in brackets is small like 1/6 and in the leading 
approximation over 1/6 one can also decouple the average of two G-s in the first term. 
Thus 

GW-^iyfu-E^GWGw) . (7) 

Our derivation of this formula in fact repeats the calculation of Green function for 
full N x N random matrices, which was done many years ago || || (see also in ]7| 
the same proof for pure band matrices). 

Due to (0) only closed chains contribute to each (V n ) in @. In terms of dual 
Feynman graphs each {V n )u may be thought as n-vertex polygon, while the averaging 
via Whick contractions may be considered as a self gluing of the edges of this polygon. 
Within this language in ((?]) we have summed up exactly all the planar (spherical) 
graphs. The nonplanar contributions lead to corrections to G of the order of ~ 1/6 



2 



for real symmetric matrices and ~ 1/b 2 for Hermitean matrices M. Nevertheless the 
equation (]?]) is exact in 1/d (||). 

It is easy to expand the equation (|6|) in the series in 1/d. Up to 1/d 2 

G{i) = - l — . (l + lU(z)G(2) 2 + G{i){AG' + £?G" 
-C/ — h[i) v 

Here prime means the derivative with respect to X = i/b and 



F(i,k) fk-i\ 1 . l^F(i,k) fk-i 
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In the leading approximation the solution of (|8]) reads: 

G°(<) = (fl - - y/(E-h(i))*-4W{ij) . (10) 

The imaginary part of G° at E = E — iO reproduces the usual in matrix models 
semicircle density. In general the 1/d 2 corrections to (|10"D may be found from (§) 
by iteration. Due to singularity of G° this corrections, which are proportional to G' 



and G", spreads to infinity when E = h(i) ± 2JW(i). Thus the nonperturbative 
treatment of equations (|) or (H) is necessary near the edge point. 

The use of expansion (H) assumes that both functions h(i) and W{%) are equally 
slow (Q). The smoothness of W is necessary because we want to consider the almost 
banded matrices, but one can sufficiently weaken the restriction for h(i). In this case 
the solution of equation fl7|) in the Breit-Wigner form is easy to found 



1 



dh 



1 

Vb 



dj 

Here the upper bound for dh/dj ensures the smallness of the regular energy intervals 
h(j + 1) — h(j) compared to the hopping matrix elements Vij (0). 

3. Turning back to the slow h(i) case let us introduce the new "scaling" variables 

E-h(0) = 2 V /W^) + ^, 

X = i/b = d x ' h x , G = -=L= + ^ , (12) 

W(0) d,b 



W(X) = W(0)+X^ = W(0) + ^ 

dX d 

h(X) = h(0)+a— , a ~ 1 . 
a 
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Here we have also made the explicit definition of the large parameter d. Substitution 
of (H) into (|) leads to 



By" + W^ 2 y 2 = e - I +a \x . (13) 



Here W = W(0). Due to (§) the term ~ AG' in (§) leads to negligible correction of 
the order of ~ 1/d 4 ^ 5 . As well the higher derivatives of G, which we have neglected 
in (HD, contribute like some powers of 1/d. 

The equation (|l"3"D is the famous Painleve I equation. Parameters B, W and a 
which still alive are of the order of 1 and may be dropped out by trivial rescaling. 
Being a part of Green the function y = y e (x), have to reproduce its analytical features. 
Thus the only singularity of y £ (x) is the cut along the axis Ime = 0. We have defined 
t ( |i~2|) as a real variable. On the other hand it is natural to consider y £ (x) as an 
analytical function of one complex variable x' = x — Now the function y(x') 

should not have singularities in the whole upper half plane. 

It is amazing that our model is not the first model of random matrices whose 
scaling limit is described by the Painleve equation. A few years ago the random 
matrix models helped to solve exactly the 2c? quantum gravity [[J. It was shown 
that the second derivative of partition function of 2d gravity with respect to the 
cosmological constant is the solution of Painleve equation. The unique solution of 
this equation, which is realized in Euclidean 2<i-gravity was found in Each 



solution of eq. (|i~3"D has an infinite set of second order poles on the complex plane 
y ~ (x — Xi)~ 2 and no other singularities. As it was pointed out by F.David fllOl 
(and as was known for mathematicians for many years) only one solution may have 
no poles in the whole half plane. This is the so called "triply truncated solution", 
which may have only a finite number of poles in the sector < Arg(x) < ^ (and 
still has an infinite set of poles within the rest ^-). Just this solution is realized in 
2o!-gravity if the matrix model is regularized via analytical continuation. Moreover 
only this unique solution can reproduce the analytical features of the Green function 
in our model @. 

The most informative quantity, which may be found from the averaged Green 
function is the single particle density 



1 VW 
p{i, E) = -ImG(i, E - i0) ~ Im y{x - - — e) . (14) 

7T 2 

It is easily seen that for real e all real solutions of (O) are unstable (have the poles 



at real x axis). The complex solutions, in accordance with (0), have all the same 



asymptotics Imy(x — > ±oo) = — ^ e — + ax). Nevertheless the nonperturbative 
imaginary part of y at negative x allows to distinguish the triply truncated solution 
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among the others [TO, ITT 



Imy(x --co) = ^ 



1 + 



exp 




'^(l + av^) 1 / 4 !— £ -x] \ ■ ( 15 ) 



Thus we see that close to the edge the single particle density behaves like 

P(hE) = -Lf(d 4/5 (E-E )) , (16) 



where / is a smooth imaginary part of the triply truncated solution. This result 
accounts exactly for the series of most singular corrections over l/d n (|],[T2|). 



Still we have not considered the corrections of the order of 1/6. These corrections 
accounts for the finite width of the band (0) and also can smooth out the singularity 
of the zero order result (fLOf). The accurate description of 1/6 effects is beyond of the 
main subject of this paper. Here we give only the result and postpone discussion for 
a separate publication |12| . For the pure band-matrix (i.e. if F(i,j) = F(\i — j\) and 



h(i) is negligible small) the density of states close to end-point reads 

p(E) = -± E <f ) (b^(E-E )) (17) 



with some smooth 0. So we can conclude, that the scaling (pT^JT3| ) should not be 
disturbed by the 1/6 corrections until d <C 6. 



4. As we have seen ( p~2] , p~3|) the variation of band width W{i) and shifting of the 
center of zone h(i) play the same role at the edge. Therefore in the following two 
sections we suppose that h = 0. 



Up to now we have supposed (|12D that the edge of spectrum E e = h(i) + 2y W(i) 
is a smooth monotonic function of i. On the other hand it seems to be interesting to 
investigate also the statistical properties of the model (@) near the bulge of E e . The 
set of scaling variables suitable for this case reads 

W(X) = W(0) -?L,X = i/b= §^d^x , (18) 



^5/6 rf 4/3 ' " ^ /W[&} ^ 7/6 dV3 



Here again W = W(0) and B are of the order of 1. Also here we give a new definition 
of large d. Simply substituting (|TjJ) into the same equation (|8|) one gets 

y" + y 2 = £ + x 2 . (19) 
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Like in monotonic case ( |i3|) the first derivative y' and all higher derivatives are sup- 
pressed by some powers of 1/d. The asymptotic solution of this equation reads 



y = -Ve + x 2 + ... (20) 
in accordance with flIQ|). For large negative e in vicinity of e + x 2 = the equation 



(19) reduces to the Painleve I and the singularity of (20) is smoothed out via its triply 



truncated solution. Thus the only new problem will be to find y at e <~ 1. 



Generally speaking the equation (|19D is not of the Painleve type. Its solutions 
also has poles on the complex plane, but each pole turns out to be also the edge of 
logarithmic cut. We are interested in the set of real x solutions y £ (x) analytical in 
the whole lower half plane of complex parameter e. Let us write down explicitly the 
real and imaginary parts of (|T9|) 




(21) 

For real energy (A = 0) the latter equation may be thought as the Schrodinger 
equation with v being the wave function, (—u) - the potential and eigenvalue equal to 
zero. Moreover it is immediately seen from (|2~0D that this Schrodinger equation could 
not have any localized solution at sufficiently large positive e. Thus the imaginary 
part of Green function (or the single particle density) vanishes in our model at large 
e. 

Our scaling equation (|19|) accounts only for the most singular corrections in each 
order over 1/d. Nevertheless, one can see, that for the exact equation (|6]) the imaginary 
part of G(i) also vanishes at sufficiently large E. In fact at large E the imaginary 
part at is small and thus it should satisfy some linear homogeneous matrix equation. 
But in general (except for the degenerate case) the linear equation has no nontrivial 
solutions. 

The exact solution of equation (|19D (or (|2l|) ) may be found only numerically. Nev- 
ertheless one can further clarify, how the imaginary part appears for real e. Suppose 
that y = Uq(x) is the real solution of ([19]) with smallest real Eq. Then the corre- 
sponding Schrodinger equation at u = uq also should have the localized zero mode 

j 4 + u 2 = e + x 2 

{ (-£-2 M0 )^ = . W 
Now it is rather simple exercise to find the solution in vicinity of Sq\ 



£o + A, y= Uo -^|^o + 0(A). (23) 
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So we see that the square root singularity of the zero order result fliCf ) which seemed 



to be smoothed out by the solution of differential equation ( |i"9|) still alive at the 
very edge of the spectrum. This result may have even experimentally measurable 
consequences. One can introduce the density of eigenstates (see (PI)) 

p(£) = $>(*,£) . (24) 

i 

This quantity is well defined because we are working at the point of maximum thick- 
ness of the wire (band). At large negative e the asymptotics of p(E) may be found 
from © 

- s > 1 , p{e) ~ (s - ef 2 . (25) 
On the other hand close to the edge 

< e - e < 1 , p(e) ~ y/e - i • (26) 

Thus the " smooth" exact result appears to be even more singular than the zero order 
approximation ( fZ5|) . 

Of course the vanishing of Jmy at e > Eq is the artifact of our approximation. 
Again in this section we have neglected all the 1/6 corrections. Comparing (|i"8|) and 
(|T7|) one can find that our results are valid if d <C b 3 ^ 5 . 

5. Up to now we have considered the averaging of only the single Green func- 
tion. The much more puzzling quantities are the correlation functions of two Green 
functions. In particular the correlators are sensitive to the effect of localization. For 
usual band-matrices the localization length turns out to be of the order of ~ b 2 |TJ 
and can not be found without nonperturbative treatment of ~ 1/6 effects. Neverthe- 
less near the thickening of the wire some sort of compulsory localization takes 
place and the eigen-functions are localized due to the "geometry" of band (wire). In 
this section we would like to find the averaged value of the squared absolute value 
of Gy at the thickening point. In fact what we can easily calculate is the so called 
"smoothed" correlator [/]], valid only for sufficiently large ImE (at least as compared 
to the interval between the individual levels). 

Like in (0) one can expand correlator in the series 



w-ii£(K(a" • 

We are mostly interested in the case of small imaginary part of E. The crucial 
observation is that in the planar limit all the Whick contractions between (V) n and 
(V) m may be done explicitly, while the sum over (V/E) s which are contracted within 
one of the multipliers in fl2"7j ) may be resumed back to G(i) or G(i)* (pi): 

> F (hh) ,^,12 F (i r 



l^| 2 = £ E \G\r-^\G\l...-^\G\) = Y,Tn{i,j) , (2f 

n=0«i,...,i„ ! ' n=0 
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where \G\ 2 = |G(i)| 2 . Like the equation (|7|) this formula accounts exactly for all ~ 1/d 
corrections. All the ~ 1/b corrections to \Gij\ 2 have been neglected in (^). Unfor- 
tunately the price for this approximation is rather high. The ~ 1/b corrections will 
blow up the correlation function at very small ImE, thus preventing the observation 
of the correlations in neighbor levels. 

The recursion formula for T-s (compare with([7|,|8])) reads 

IWU) = J2\G\l^j^T n (k,j)= (29) 
k 

= W\G\ 2 T n + A\G\ 2 T' n + B\G\ 2 T^ + ... , 

where A, B, \G\ 2 are taken at the point i and the prime means the derivative with 
respect to X = i/b. Now if one uses the scaling variables (|T^) and also introduces the 
additional variable (time) instead of n the differential equation for T(t,x,Xj) is easy 
to write down 

n 

t= d^> T = T" + 2ReyT . (30) 

Here and below we choose W(0) = B = 1 for simplicity. It is seen immediately 
from fl2~T|) that T spreads to infinity if Im e — > due to zero mode of the Schrodinger 
equation in the r.h.s. of (|30"D. 

For small A = —line it is convenient to consider the linearized version of the 
system (|2~lf) 

u = Wo + / , v = v + g , /, g ~ A , 
/ f" + 2u f = 2v g 

\ g" + 2u g = -A -2v f . [6i) 

Here u , v are the solutions of (pl|) for A = and hence v is also the zero mode 
of operator ^-y + 2mq. Thus one can write down the consistency conditions for two 



equations 



2j v 2 g = , X J v = -2 J v 2 a f . (32) 
The solution of (|30|) may be found as a sum over eigenmodes of the operator: 

d 2 d 2 

- — -2u = -_-2m -2/ , 

ctar dx z 

r(t,x,o;j) = fcjTJsYl \n)(n\exp(-e n t) , (33) 

where e n and |n) are the eigenvalues and normalized eigenmodes and the overall factor 
in T was found from the initial condition (see fl28|) ). For small A and large t only the 
almost zero mode (n = 0) survives in (|33|), which is: 

|0) -^' e °- _2 T^- A M- (34) 



s 



Finally the "smoothed" pj correlator reads 



]7t-[2 _ 1 Vo(xi)v (xj) ^r r r^_ g(vg) 

Here we use that A = d^ 3 ImE (]T8|) . We have also shown the sum rule, which this 
correlator evidently satisfy. Function vq(x) may be found numerically as a solution 
of the simple system of differential equations (PI). 
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